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Abstract. We present a study of the temperature (T) dependence of the dc electrical 
conductivity of polyaniline pellets doped over a wide range. A crossover between low- 
T variable-range hopping and high-T nearest-neighbor hopping has been found below 
room temperature for most of the samples, which is atypical for conducting polymers. 
This allows us to construct an experimental diagram that accounts for different regimes 
in the hopping electron transport and closely resembles predictions of the Foglcr-Tcber- 
Shklovskii model in three dimensions. 
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Particle transport through a disordered potential consisting of deep, localizing 
potential wells has for decades been puzzling. This transport usually occurs via 
thermally assisted hopping from one localized site to another, good examples being 
electron transport in strongly disordered conductors [T], and vortex creep in type-II 
superconductors [2]. If the well depth varies in space strongly, a particle hops to a 
site determined by the most favorable balance between the hop length and the energy 
difference between the localized states involved. When thermal activation is strong 
enough to diminish the importance of the spatially varying well depth, the above 
mechanism - called variable-range hopping (VRH) - is replaced by nearest-neighbor 
hopping (NNH). The most suitable observable for investigating these phenomena is 
electrical conductivity, which is a linear-response quantity that can contain enough 
information on the nature of a hopping electron transport (HET). This approach has 
been dominant ever since the early works of Anderson [3] and Mott [I] on localized 
electron states in solids and their thermally assisted derealization. Mott argued [I] 
that HET electrical conductivity a generally depends on T as 

a(T)= Va ex V [-(T a /T) a }, (1) 

where rj a ,T a , and a < 1 contain information on the underlying microscopic mechanisms. 
These parameters may depend on dimensionality (d), interactions relevant for the 
transport, disorder level, localization length of electrons, etc. For VRH, a < 1, whereas 
NNH obeys the Arrhenius law (a = 1). 

VRH of noninteracting electrons results in the Mott law a = (1 + d)~ . Coulomb 
interaction leads to an \E — Ep\^ energy (E) dependence of the electron density of 
states (DOS) close to the Fermi energy Ep, with fi being a constant. This opens a 
soft (Coulomb) gap in the DOS, and results in the Efros-Shklovskii (ES) law where 
a = 1/2 is universal [6]. The above models address isotropic systems, e.g., disordered 
semiconductors. However, HET occurs in other materials as well, of which we pay 
attention to conducting polymers. For many of them, Eq. ([T]) accounts for experimental 
data remarkably well [7], with a = 1/4 and a = 1/2 appearing often. Other a < 1 have 
been found as well, including a = 2/5 [SI El EU El EE2] that cannot be explained 
neither by the traditional models (Mott, ES) nor by similar approaches [131 HH HS1 116] - 
NNH is, on the other hand, much less common in conducting polymers, and a = 1 is 
observed [UJ HH [19] infrequently, probably because conditions for its occurrence are 
seldom established at T lower than the material degradation temperature. 

A model with a potential to account for the complexity of cr{T) in conducting 
polymers was developed by Fogler, Teber and Shklovskii (FTS) [50] . They studied the 
HET in coupled chain-like conductors in the presence of disorder-dependent Coulomb 
interaction potential that affected DOS similarly as in the ES model. The FTS model 
predicts a = (1 + fi) / (1 + // + d), which for strongly coupled chains (d = 3) and 
fi = 0, 1, 2 leads to a = 1/4, 2/5, 1/2, respectively. At low T, highly disordered samples 
should exhibit a = 1/2, and more ordered samples a = 2/5. In the former case, the 
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transition to NNH as T increases should be direct, i.e., a — 1/2 — > 1, whereas in the 
latter case, one should observe a = 2/5 — > 1/4 — > 1. On the basis of these considerations, 
FTS proposed a d — 3 diagram where crossover temperatures between different a values 
were plotted vs suitably parameterized disorder level. A support to the FTS predictions 
has come from recent experiments on VRH in polyaniline (PANI) [HI HZ] which is a 
well-known polymer that becomes conducting when doped via protonation by an acid 
In HCl-doped samples (PANI-HC1), a exhibited crossovers between 1/4, 2/5, 
and 1/2 [TTJ, but the transition from VRH to NNH below room temperature (RT) was 
identified only occasionally [17]. Hence, the FTS diagram - which is in many respects 
central to the model - could not have been constructed from the experimental data. 
A similar VRH behavior has been found in samples of PANI-DBSA, i.e., PANI doped 
with dodecylbenzenesulfonic acid (DBS A) [12]. In this case, however, VRH-to-NNH 
transition occurs well below RT over a wide doping range, and this opens a possibility 
to explore the viability of the FTS diagram. 

In this paper, we present a systematic study of the VRH and NNH in pressed 
PANI-DBSA pellets doped over a wide range. For the first time, we construct the FTS 
diagram from experimental data and confirm its correctness. Data for those of our 
PANI-HC1 samples which show NNH below RT supplement the diagram nicely. Hence, 
our findings may have consequences for general viewpoint on the HET in conducting 
polymers and similar conductors. 

Detailed information on the sample synthesis, doping, and characterization 
(nuclear spectroscopy for elemental analysis, X-ray diffraction for structure, magnetic 
susceptibility) can be found in [23], and here we list the main conclusions. Q = 4ns/ W, 
where ns and are concentrations of sulfur and nitrogen in PANI-DBSA, accounts for 
the doping level and can be as high as 3.39 (while Q = 2 corresponds to the complete 
protonation). This means that the amount of DBSA in a sample can exceed that of 
the protonation sites on PANI molecules. Measured transport and magnetic properties 
imply that all DBSA participates in the protonation when Q < 2, whereas Q > 2 means 
that PANI is completely protonated and extra DBSA is captured between the protonated 
PANI molecules. For all Q, the error of which is at the level of a few percent, no signature 
of crystalline ordering is found, a of the samples increases monotonically with increasing 
Q over the whole doping range. While a (Q) for Q < 2 is a straightforward consequence 
of the protonation, the increase at Q > 2 has been attributed to a rearrangement of the 
fully protonated PANI molecules into a configuration where interchain electron hops 
are enhanced [12]. Dc conductivity measurements were carried out in a closed-cycle 
refrigerator from the lowest achievable temperature of ~ 10 K up to RT. 

Experimental er(T) curves for the obtained range of Q are in Fig. Q](a) plotted by 
the symbols as log cr vs T~ 1//2 , which provides good insight into low-T data. Below the 
dashed line, the plots are obviously straight. Above the dashed line, the symbols and 
the corresponding straight solid lines (guides to the eye) do not match well, hence there 
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Figure 1. (a) Experimental log a vs T -1 / 2 (symbols). The solid lines are 
[here and in (b)] guides to the eye. (b) Data above the dashed line in (a), plotted 
as log a vs T~ 2 ' 5 (symbols). Inset to (b): logtr vs T- 1 / 4 for Q = 3.39, exhibiting 
linearity for 110 < T < 145 K. The corresponding a determined numerically 
[as exemplified in (c)] is 0.27. (c) log VP vs logT of a selected sample (symbols), 
displaying two distinct ranges of linearity with a crossover at T*. The solid lines 
represent linear fits that yield a numerically, a resulting from this procedure is listed 
in the legends to (a) and (b). 
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is no linear dependence. These data are in Fig. QJb) plotted by the symbols as log o vs 
T -2 / 5 , which results in linear plots. From Figs. Ufa) and[IJb), we can conclude that Eq. 
(JTJ) accounts for low-T data adequately, with a = 1/2 for low Q, and a = 2/5 for high 
Q. A similar behavior was found for PANI-HC1 (in which Q = 2 cannot be exceeded |j) 
but the change from a = l/2toa = 2/5 occurred at a Q lower than here [TTj . 

It is difficult to note in Figs. GJa) and QJb) that the above low-T values of a do 
not hold for the high-T data in narrow ranges on the left-hand side of the plots. As we 
show later, the proper exponent in this case is a = 1 irrespective of the low-T value, 
indicating that NNH is the mechanism that underlies the HET in this regime. Moreover, 
for the sample with the highest Q = 3.39, there is a narrow T range (~ 110 — 145 K) 
where a = 1/4, as shown in the inset to Fig. QJb). Therefore, as T increases, a of this 
sample changes as 2/5 — )■ 1/4 — >• 1, in contrast to the other samples where this change 
is either 1/2 — > 1 or 2/5 — > 1. 

Although logo" vs T~ a plots are indicative with regard to a, we carried out an 
alternative analysis as well. Application of W = d In a/d In T to Eq. ([T]) leads to —a 
being the slope in the resulting linear log W vs log T dependence, which enables a more 
precise extraction of a from experimental cr(T). This is demonstrated in Fig. QJc) for 
a selected sample. Two distinct ranges of linearity are discernible, corresponding to 
VRH (a = 0.42) and NNH (a = 1.01), as indicated (for all our data, error in thusly 
determined a is < 10%). T* is a crossover temperature between the VRH and NNH. 
Crossover at T* is generally not sharp, and the T range between regimes of well defined 
low-T and high-T values of a can be up to ~ 40 K wide. A large crossover width can 
be an obstacle for identifying a = 1/4, because this VRH behavior should - according 
to the FTS model - take place only in a narrow T range below T*. Our result for 
the Q = 3.39 sample, where a = 1/4 appears only over 35 K below T*, supports this 
prediction. Moreover, we cannot rule out that a = 1/4 in some other samples might be 
hidden due to a large width of the transition at T*. 

In Fig. [21 we use loger vs T" 1 plots to show the high-T part of cr(T) at Q > 0.38, 
for which T* is lower than RT. Larger symbols mark T* which decreases with increasing 
Q. We also note that magnetic susceptibility x of the samples undergoes a significant 
change [jj, similar to that described in [23] for PANI-HC1, at basically the same T*(Q) as 
here. The fact that T* leaves a signature both in a and x implies a qualitative change 
in the electronic system at T*. 

From the above analysis, it is rather clear that several predictions of the FTS 
model agree with our experimental results. The remaining issue is whether we can 
construct the FTS diagram - the proposed shape of which is sketched in the inset to 

|For PANI-HC1, Q is calculated by replacing ng with CI concentration tiq. In fully dry samples, 
no extra HC1 can be intercalated between PANI chains. 

§M. Bacani. I. Kokanovic, M. Novak, and D. Babic, unpublished 
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Figure 2. Experimental log a vs T _1 (symbols) at Q > 0.38, for which T* is below 
RT. The solid lines are guides to the eye. T* for each sample is marked by a larger 
symbol. Numerical values of a extracted numerically are listed in the legends. 



Fig. [3] - from our data. This diagram maps areas of different a in a T versus disorder 
level plane. VRH exponents 1/4, 2/5, and 1/2 can all appear only for d = 3 [20] . 
and this dimensionality should hence apply to our samples. Since data on crossovers 
between different a are available from our cr(T), constructing the FTS diagram reduces 
to parameterizing the disorder level. This is, however, not straightforward. Disorder 
in conducting polymers arises from a combination of molecular-scale disorder (see e.g. 
[25] and references therein) and structural inhomogeneities at mesoscopic scale lengths 
[26J. Thus, information on averaged structural ordering (e.g., from X-ray diffraction) 
can be instructive but not conclusive in inferring what level of disorder a hopping 
electron experiences, particularly in the case of our PANI-DBSA samples which are 
fully amorphous. Having this in mind, we approach the problem by using a macroscopic 
property of cr as a parameter accounting for disorder. 

We believe that rj a for a = 1, i.e., experimental rji, can be used as an appropriate 
measure of the disorder level in our samples. This plausible assumption, as we show later, 
leads to our results reproducing the FTS diagram, and the reasoning behind choosing 
rji is as follows. Generally, the exponential part in Eq. (P) pertains mainly to thermal 
aspects of a while r\ a is related to non-thermal effects - in which disorder plays a crucial 
role. Since a oc rj a , as disorder level decreases, rj a should increase. This holds although 
rj a also contains a contribution of the charge injected by the doping. Protonation 



Hopping Electron Transport in Doped Polyaniline: an Experimental Verification of the FTS Model 



7 



not only increases the charge- carrier density (which also affects the screening of the 
Coulomb interaction and thusly influences T a [IT]) but also reduces disorder in the 
hopping landscape [27]. This occurs because the Coulomb interaction energy between 
dopant anions - which are bound to the protonated sites together with the potential wells 
- is lowered in a more ordered configuration. Another possible effect of adding a dopant 
is enhancement of the interchain hopping pj2] , which is also a process that homogenizes 
the potential. Of all r] a , we focus on r\\ because it accounts for the described processes 
in the same way for all the samples which exhibit NNH, regardless of the low-T value (s) 
of a. 
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Figure 3. Areas of different values of a, shown in a T vs 771 plane. The open 
symbols correspond to the PANI-DBSA samples, and the solid ones to two PANI-HC1 
samples that exhibit a — 1 below RT. The lines are guides to the eye. The uppermost 
symbols correspond to T* , whereas the bottom symbols show Tlow which is the lowest 
temperature at which a was measured. The intermediate symbol for the two highest 
771 values mark crossover between a — 2/5 and a = 1/4. A sketch of the theoretically 
predicted diagram [20 is shown in the inset. 



In Fig. El we plot areas of different values of a in a (t)i,T) plane, as indicated. 
Each rji value is related to one sample, and the symbols have been extracted from the 
corresponding cr(T). The vertical distance between two symbols for a given rji is the 
T range of the appearance of the indicated VRH a. The open symbols correspond 
to the PANI-DBSA samples, and the solid symbols to the only two of our PANI-HC1 
samples (both with Q < 2) that exhibit well defined T ranges of a = 1 below RT 
(their undoped PANI being the same as for the PANI-DBSA samples). The uppermost 
symbols correspond to T* (the error bars indicate the crossover widths), whereas the 
bottom symbols represent the lowest temperature Tlow at which a was measured (the 
low-T limit of the cryogenic system, or T at which a is immeasurably small). For the 
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two samples with the highest rji, the intermediate symbol represents crossover between 
a = 2/5 and a = 1/4. The dashed lines are guides to the eye, depicting characteristic 
temperatures extracted from cx(T), whereas the dotted lines mimic other lines from the 
theoretical diagram. 

The similarity of the experimental and theoretically predicted diagrams is obvious. 
Furthermore, the two PANI-HC1 samples show the same behavior as the PANI-DBSA 
samples. We believe that this is not a coincidence, because the starting undoped PANI 
has been the same for both classes of samples and every processing (set by the choice 
of the dopant and the doping method) modified the same starting material. On the 
other hand, a particular processing leaves a distinct signature in rji. For instance, 
both PANI-HC1 samples are almost fully protonated but they were processed differently 
and hence do not display the same properties [T7]. The sample with rji ~ 2000 
S/m exhibits a = 1/4 while this regime is absent for the sample with wi ~ 1100 
S/m. Moreover, the PANI-HC1 samples are partly crystalline [UJ while the PANI- 
DBSA samples are amorphous, which additionally supports the conclusion that degree 
of crystallinity is not a proper measure of the (dis)order in a potential that hopping 
electrons experience. It seems, however, that the doping-induced ordering of the hopping 
potential must be strong enough for different classes of doped PANI to show the same 
HET behavior quantitatively, since our less ordered PANI-HC1 samples from the same 
starting material do not exhibit NNH below RT [11]. That is, the wi < 1000 S/m parts 
of the FTS diagrams for PANI-HC1 and PANI-DBSA differ quantitatively although the 
corresponding VRH properties are similar [11]. In contrast to PANI-HC1, T* for PANI- 
DBSA does not increase above RT until rji decreases below ~ 75 S/m. This increase is 
depicted by the dotted line for rji — > 0, which can be drawn thusly because there is no 
NNH below RT for the least conducting of the PANI-DBSA samples (see Fig. [1]). 

In developing their model, FTS primarily had in mind systems such as stripe phases, 
quantum wire arrays, carbon-nanotube films, compounds like Bechgaard salts and Fabre 
salts, etc. They also mentioned [281 EH] conducting polymers as candidates since VRH 
has been reported extensively for these materials. Yet, FTS were seemingly not fully 
confident in the applicability of their model in this case, possibly because of the intricate 
underlying morphology. It turns out, however, that the model is robust enough to 
account for these complex HET conductors as well. 

In conclusion, the cr(T) curves (between 10 K and RT) of our doped PANI pellets 
reveal both VRH and NNH. This is observed for doping with DBSA over a wide range 
as well as for nearly full protonation by HC1. The two HET mechanisms extend over 
T ranges that are sufficiently large to permit studying of the full complexity of the 
FTS diagram. We determine crossover temperatures between regimes of different HET 
exponents. When we plot these temperatures against a measure of the disorder level, 
that is, prefactor of the NNH Arrhenius plots, we obtain a remarkable agreement with 
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the predicted shape of the FTS diagram. This implies that the FTS model is both 
correct and applicable to conducting polymers. 
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